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ABSTRACT 

We present kinematic simulations of a galactic dynamo model based on the large 
scale differential rotation and the small scale helical fluctuations due to supernova ex- 
plosions. We report for the first time direct numerical simulations of the full galactic 
dynamo using an unparameterized global approach. We argue that the scale of helicity 
injection is large enough to be directly resolved rather than parameterized. While the 
actual superbubble characteristics can only be approached, we show that numerical 
simulations yield magnetic structures which are close both to the observations and to 
the previous parameterized mean field models. In particular, the quadrupolar symme- 
try and the spiraling properties of the field are reproduced. Moreover, our simulations 
show that the presence of a vertical inflow plays an essential role to increase the mag- 
netic growth rate. This observation could indicate an important role of the downward 
flow (possibly linked with galactic fountains) in sustaining galactic magnetic fields. 
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1 INTRODUCTION 

It is widely accepted that magnetic fields of planets, stars 
and galaxies are generated by dynamo action, i.e., by the 
magnetic field amplification due to electromagnetic induc- 
\ tion asso ciated to the motion of an electrically conduct- 
ing fluid (|Moffattlll978h . The flow of gas in the interstellar 
medium appears to convey the essential ingredients for suc h 
dynamo action (|Wielebinski| [T990: Ros ner fe Demcalll989h . 
Differential rotation in the galactic disk creates a strong 
shear along the radial direction. This shear is very efficient at 
stretching radial magnetic field lines in the azimuthal direc- 
tion (this is known as the o;-effect). In combination with this 
large scale effect, the turbulent motions at small scales pro- 
vide a cyclonic flow generating poloidal magnetic field (this 
is the so-called a-effect). Together, both effects suggest the 
possibility of an a — u) type of dynamo that mig ht be respon- 
sible for generating the galacti c magnetic field (|Parkerlll97ll ; 
IVainshtein &; Ruzmaikinlll97lh . Let us note that alternative 
models for dynamo action in galaxi es have been prop osed 
through the action of cos mic rays Jl lanasz et al.1 2004) or 
in a cosmological context ( Wang &; Abell [2Q07h , which will 
not be discussed here. The apparent scale separation be- 
tween the shear and the turbulent motions has often been 
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invoked to introduce a mean field approach for the galactic 

i II 1 | ; it 1 1 i 

dynamo (Beck et al. 1996; Ferriere 1998). In such formalism, 
an equation for the large scale magnetic field only is solved, 
the effect of small scales being parameterized by an a term 
([Krause &; Raedlerlll980h . Relying on mean field equations 
has proven to b e a very efficien t approach to the galactic dy- 
namo problem (Ferriere 1992). It is for example an efficient 
way to achieve moderate simulation time. However, the re- 
sults of mean field simulations are intrinsically limited by 
strong assumptions such as scale separation or the statisti- 
cal properties of turbulence. It is thus interesting to study 
galactic dynamos with direct simulations of the full problem 
by properly treating the small scale flow associated with the 
turbulence in the interstellar medium and thus solve for the 
magnetic field at all scales. 

It is often assumed that the most importance source 
of turbulence in t he interstellar medium comes from super- 
nova explosions ([McCrav &; Snow| [l979). The positions of 
these explosions are not completely random in the disk but 
they often occur in cluster. This produces giant expanding 
cavities of gas known as superbubbles. These explosions oc- 
curring in a rotating galaxy, the expansion is affected by 
a Coriolis force. This yields cy clonic motions and thus a 
strong helicity in the gas flow (|Ferrierdll998h . In such a 
framework, however, it is worth noting that the scale sepa- 
ration mentioned above is not dramatic. Superbubbles have 
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Figure 1. Field lines of the magnetic field (blue) are represented for Rm = 7000. Both the spiral structure of the magnetic field and its 
quadrupolar symmetry can be identified. Isosurface corresponding to 1% of the peak kinetic energy are also represented (red). 



typi cal sizes of the ord er of a few hundreds parsecs (pc, 
see lOev fe Clarkel Il99/t ). This is smaller, but not dramat- 
ically smaller than the typical vertical scale of the galaxy 
(rsj 1 kpc). Given modern day computational resources, these 
numbers suggest that direct numerical simulations (i.e. nu- 
merical simulations that do not rely on an ad hoc parameter- 
ization of the small scales, for example through the a-effect) 
are within reach. Indeed, (|Gressel et al.l l2008) recently pre- 
sented such simulations. To cope with the large resolution 
still needed to address this problem, they adopted a local 
approach based on the shearing box model. Their results in- 
dicate a good agreement of the local approach with mean 
field models. However, the local approach they used pre- 
cludes any global diagnostics, such as the global structure 
of the field, to be established. 

The purpose of this letter is to present such global nu- 
merical simulations, resolving the magnetic field at all rel- 
evant scales in the galaxy (i.e. from 100 parsecs to 10 kilo- 
parsecs). To reduce the computational burden that would 
be associated with full MHD simulations, we work in the 
kinematic regime: we solve the induction equation using a 
prescribed and time dependent gas flow. The later is set by 
using an analytical velocity field which intends to reproduce 
the large scale shear associated with rotation and the effect 
of superbubble explosions on the interstellar medium. 



2 NUMERICAL MODEL 

The direct numerical simulations presented in this letter 
are fully three dimensional. We solve the induction equa- 
tion governing the evolution of the solenoidal magnetic field 
B in a cylindrical coordinate system (r,<fr,z): 

<9B 

— = V x (u x B) + Rm- 1 AB , (1) 



written in dimensionless form using the advective timescale. 
The magnetic Reynolds number Rm is defined as Rm = 
fioaroUo, where /io is the permeability of vacuum, the typ- 
ical length scale ro = 10 kpc is the radius of the galactic 
disk, the typical velocity scale Uo is the velocity of the large 
scale flow (i.e. the differential rotation of the galaxy) and 
a is the conductivity of the plasma 0. In our simulations 
the vertical extent of the galactic disk is H = ro/10, z thus 
range from — 0.05ro to +0.05ro- We restrict our attention 
here to the kinematic problem, ignoring the back reaction 
of the magnetic field on the flow. The velocity field u used in 
Eq. (pQ) is analytical and represents the differential rotation 
of the galaxy and the supernova explosions. This approach 
also means that we do not explicitly consider important ef- 
fects such as density stratification in the vertical direction 
or the induction effect which wou ld be due to interstellar 
turbulence (|Ruzmaikin et al.lll988l ). 

Eq. (pQ) is solved using a finit e volume approach. The 
method is described in details by ([Tevssier et al.l [20061 ): it 
uses the MUSCL-Hancok upwind method. The solenoidal 
character of the magnetic field B is maintained through the 
const rained transport algorithm (Yee 19661; [Evans &; Hawlevl 
1 19881 ). We rely here on the so-called pseudo- vacuum bound- 
ary conditions for the magnetic field. This corresponds to 
imposing B x n = at all boundaries of the computational 
domain. These boundary conditions are not fully realistic, 
but they are often used in parameterized models of galactic 
dynamos and simple to implement. These boundary con- 
ditions are known to modify quantitative results (such as 

1 We could have adopted alternative definitions of the Reynolds 
number, for example Rm' = fj,oaHV s /2, where V s is the sound ve- 
locity, which is equal to the terminal velocity of the superbubbles 
(see later in the text). With this definition, Rm' = Rm/200, and 
the maximum value achieved in this work would be Rm' = 500. 
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the threshold value for dynamo action) b ut not the global 
qualitative solution (|Gissinger et al.ll2008h . We now turn to 
a detailed description of the velocity field being used. It 
is the sum of two terms: rotation around the vertical axis 
and modification of the flow by superbubbles. In our sim- 
ulations, we use the following prescription for the rotation: 
U = Uo e$ , with a constant Uo . This is a good approxi- 
mation since the angular velocity is observed to be roughly 
proportional to 1/r in galaxies. The effect of supernova ex- 
plosions is more subtle to implement. We decided to con- 
sider the effect of superbubbles only and ignore here iso- 
lated supe rnovae, as the e nergy input of the former is largely 
dominant (|Ferrierelll998l ). Considering superbubbles rather 
than smaller isolated supernovae yields larger scales which 
directly translates into resolutions affordable with modern 
days computing resources. Let us consider first the explosion 
of one superbubble, in a local sphe rical co ordinate system 
(r',9' ,</>'). Following the work of (jFerrierd \19m ) . we work 
under the simplifying assumption that each explosion rem- 
nant has a perfectly s pherical shape. We thus use the simple 
radius evolution law ([Weaver et al.|[l977h : 



At" 



(2) 



During the expansion of each superbubble, the rotation of 
the galaxy yields a Coriolis force which tends to deflect the 
initially radial expansion and create cyclonic motions. This 
is an essential step in cl assical mean field a — u) description 
of the galactic dynamo (|Ferrierelll998h . This Coriolis effect 
can be evaluated by solving the equation of gas motion: 



2CI x v r 



(3) 
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where F e is a force leading to the radial expansion described 
by Eq. (|2]). Integrating Eq. © in the radial direction leads 
to the expansion (|2]) . The azimuthal velocity is obtained by 
integrating the equation J3} in the azimuthal direction. In 
doing so, we made the approximation that the Coriolis force 
on the superbubble is only due to the radial expansion of 
the shell. Inside the superbubble, we assume a linear vari- 
ation of velocity in radius. An important parameter is r' c , 
the critical size reached by the superbubble for which the 
pressure in the cavity becomes comparable to that of the 
surrounding medium. At this point, we consider that the 
bubble merges with the interstellar medium. This situation 
generally occurs when the radial velocity of the shell become 
comparable to the velocity of sound in the medium. In our 
modeling, this critical velocity numerically determines the 
end of existence of a superbubble. The velocity field associ- 
ated with a superbubble therefore vanishes when the radial 
velocity reach this critical velocity v c . This radial expansion 
and the associated Coriolis force totally determine the flow 
at small scales. In most observed galaxies, the spatial dis- 
tribution of explosions in the galaxy is rapidly decreasing 
away from the midplane of the disk. For simplicity, we will 
assume here that all explosions occur in the midplane only, 
but with random position in the disk. In actual galaxies, 
there is a large observed dispersion of data about superbub- 
bles, yet averaged values for the expl osions rate of superbub- 
bles are fo — 4.5. 10 -7 kpc _2 yr _1 (El megreen &; Clemens! 
1985). Such parameters, however, are still out of reach of 
present computations (especially because of the high explo- 
sion rate which implies large numbers of superbubbles to 
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Figure 2. Growth rate of the magnetic energy as a function 
of the magnetic Reynolds number. Black circles correspond to 
simulations without infall (7 = 0), whereas the filled red square 
corresponds to a simulation with 7 = 0.03. 



be handled at the same time). We use here a lower rate of 
superbubbles, but more powerful explosions, thus leading 
to a similar helicity input. In the simulations reported here, 
/ = /o/50, r' c = 0.4, A = 0.35 and v = 0.6. This corresponds 
to about 150 superbubbles expanding in the galactic disk at 
a given time in the simulations. In some cases, we will also 
take into account a downward flow. Due to the simplicity 
of our model, t his velocity could be attrib uted to turbulent 
diamagnetism (|Sok oloff & Shukurov 1990) or to the galacti c 
fountain mechanism (Shapiro & Field ll976l :lB regman 1980). 
As an attempt to describe these effects, we add the following 
vertical velocity to the flow: 



v x (z) 



-jz 



2 /2(3 2 



2tt2/3 



(4) 



It is antisymmetric with respect to the midplane and van- 
ishes for z = 0. Moreover, the infall velocity decreases far 
away from the midplane. The parameter f3 controls the ex- 
tension of the infall region and we use here f3 — r' c /?> so 
that the maximum of the infall is near the region where su- 
perbubble explosions tend to accumulate the matter. 7 is 
a free parameter controlling the amplitude of the vertical 
velocity. We will use here 7 = 0.03 throughout this paper 
corresponding to a typical velocity of 6 km.s -1 . Despite the 
simplifications implied by working in the kinematic regime, 
large spatial resolutions are still needed in order to correctly 
describe the evolution of the superbubbles at small scales. In 
the runs presented here, we used a resolution of N r = 200, 
N 6 640, Nz 36. 



3 RESULTS 

3.1 General features 

We performed seven simulations for different magnetic 
Reynolds number ranging from Rm = 100 to Rm = 10 5 
(this would correspond to Rm' between 0.5 and 500). We 
choose to stop the simulations after a few resistive times, 
when the growth rate of the magnetic energy is statistically 
invariant and the exponential growth is well established. 
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Figure 3. Structure of instantaneous magnetic field in (r, </>)- 
plane in the midplane of the galaxy at z = 0. Magnetic field lines 
projected in the (r, </>)-plane are represented by black lines and 
the color code reflects the strength of Bz. 




Figure 4. Magnetic field just before the top of the domain. Note 
that the magnetic field is smooth due to the weak effect of the 
superbubbles at this altitude. The sign of is reversed compared 
to the midplane. 

For all of these simulations, we measure the growth rate 
of the magnetic energy. It is displayed on figure [2] as a func- 
tion of the magnetic Reynolds number Rm. It is negative 
when the magnetic Reynolds number Rm is smaller than 
Rm c ~ 500. It is positive for larger Rm, indicating exponen- 
tial amplification in that case. For Rm = 10 5 , the growth 
rate is a = 0.6 Gyr -1 . Such growth rates are comparable 
to the ones obtained by (|Gressel et alj|2008l ), although they 
seem to be larger in our case. 

The result of a typical simulation (Rm — 10 5 ) once the 
exponentially growing phase is reached is illustrated in figure 
[1] which shows simultaneously the structure of the magnetic 
field and that of the flow. Many superbubbles (red isosur- 
faces) are present at a given time in the model. We also 
show field lines (plotted in blue) of the magnetic field. The 
observed magnetic structure is the results of the combined 
effects of the superbubble explosions and the differential ro- 
tation of the disk. The colored slice shows the magnetic en- 
ergy in the equatorial plane. It is strongly fluctuating due 
to the complicated nature of the flow. The overall topology 
of the magnetic field is complex. We now turn to a detailed 
study of its structure. 



3.2 Structure of the magnetic field 

The structure in the (r, 0)-plane is complicated and varies 
with the altitude z. Figure [3] shows the magnetic field in 
the midplane of the galaxy (the solid lines represent field 




Figure 5. Magnetic field at z = zO/2. Note that B^ change sign 
when the radius is increased. 
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Figure 6. Magnetic structure of (^-averaged magnetic field in the 
(r, z)-plane for: Rm = 1000 (a.), Rm = 3000 (b.) and Rm = 10 4 
(c. and d.). Flux expulsion is clearly visible on plot c. Plot d. cor- 
responds to a simulation including vertical infall. Flux expulsion 
can thus be counteracted. 



lines projected in this plane and the color code indicates the 
strength of Bz). 

In the midplane, the field is organized in a spiral struc- 
ture. The sign of B^ is constant along the radial direction. 
Near the axis of rotation of the disk, the azimuthal field 
largely dominates all others components, but rapidly goes 
to zero at the inner boundary in order to satisfy the bound- 
ary conditions. At larger radii, the vertical field is negli- 
gible while B^ and B r are now comparable. Their rela- 
tive value is given by the magnetic pitch angle, defined as 
Pb = atan(B r / B<f>). It is remarkable that, despite the fact 
that numerical parameters are far from actual values, ps is 
very close to the observations: except near the unrealistic 
boundaries of the domain, the pitch angle is in general close 
to 15°, which is in agreem ent with the r ange [—30°,— 10°] 
observed in real galaxies ()Shukurovl 120071 ) . An average of 
the pitch angle in radius from r = 0.2 to r — 0.8 gives 
Pb — —15°. This is also in agreement with (jGressel et al.l 
2008) although smaller, as they report a pitch angle around 
10°. 

At higher altitudes, the structure of the field is much 
more complicated. By increasing z, we observe that B^ can 
change sign. We always observe opposite sign of B^ between 
the midplane (z = 0) and the halo (z — +zq) (see figure 2]). 
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At intermediate altitudes, can also reverse sign along the 
radial direction itself, as it is shown in figure [5] 

While the structure in the (r, 0)-plane is not very sen- 
sitive to the resistivity, the magnetic field in the (r, z)-plane 
presents different behaviors depending on the value of the 
Reynolds number Rm, as shown in figure [6] A quadrupolar 
structure is ubiquitous in all simulations but the location of 
the magnetic loops does depend on Rm. Indeed, the effect 
of superbubbles is located near the midplane and produces 
strong expulsion of magnetic field in the halo of the galaxy. 
For weak Rm, the magnetic resistivity counteracts this effect 
through vertical diffusion. For larger Rm, the weak mag- 
netic resistivity cannot balance anymore the strong vertical 
expulsion of magnetic field due to multiple explosions. As 
a consequence, the quadrupole becomes unrealistically con- 
fined to the halo of the galaxy (Fig. [6]:), far from the active 
region. Although the creation of this external shell does not 
totally inhibit dynamo action, it clearly decreases the mag- 
netic growth rate. 

3.3 Effect of vertical infall 

This behavior indicates how the diffusion of magnetic field 
can play two opposite roles: on the one hand, it is obviously 
defavorable to dynamo action by increasing resistivity in the 
induction equation. On the other hand, diffusion can be fa- 
vorable by preventing magnetic flux expulsion away from 
the midplane region where the small scale flow is important. 
However, for weak resistivity (as is the case in real galax- 
ies) , superbubbles expel the magnetic field out of the active 
region of galactic disk, thus inhibiting dynamo action. In 
that case, adding a vertical inflow by using Eq. Q proved 
to be an essential ingredient to dynamo action. This vertical 
flow indeed pumps the magnetic field from the halo to the 
midplane, which increases considerably the growth rate of 
magnetic energy. For Rm — 10 4 for example, the growth rate 
increases from a = 0.7 Gy _1 without inflow to a = 1.4 Gy _1 
with vertical inflow (filled red square on figure [2]). As seen on 
figure [fiji, the magnetic field structure is again quadrupolar 
in that case and is spread out over the whole galaxy. 



4 CONCLUSION 

We have shown that according to our simple model, it is 
possible to perform numerical simulations of the galactic dy- 
namo without the need for a mean field formalism. We thus 
avoid assumptions in the scale separation and can control 
more rigorously the origin of the source term in the induc- 
tion equation. Our simulations yield magnetic field with two 
main characteristics: a quadrupolar symmetry in the (r, z)— 
plane and a roughly axisymmetric spiral configuration in 
the (r, (/>)-plane. Both characteristics are in good agreement 
with observations and confirm previous studies that used a 
mean field approach. A detailed study of the magnetic field 
topology shows a complicated structure, with reversals of 
along the radial or vertical directions. Another inter- 
esting features of the present work are the paradoxical role 
of superbubbles in the limit of very weak magnetic diffu- 
sion. Indeed, the turbulent flow due to explosions is, with 
the differential rotation, an essential ingredient of the a — uj 
dynamo but also inhibits dynamo action by confining the 



magnetic field in the halo of the galaxy. In this context, the 
vertical inflow of interstellar gas appears as the third main 
ingredient needed for dynamo action. The downward flow 
observed in galaxies could thus be an essential mechanism 
of galactic dynamo theory. 
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